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Abstract 

The projected random phase approximation (PRPA) for charge-exchange excitations is de- 
rived from the time-dependent variational principle. Explicit results for the unperturbed 
energies (including the self-energy corrections), the PRPA matrices, and the transition ma- 
trix elements are presented. The effect of the projection procedure on the two-neutrino (3j3 
decay in "^^Ge is briefly discussed. 



^Fellow of the CONICET from Argentina. 



In a recent work |1| simple formulae have been designed for the 0^ — > 0+ double beta 
decay matrix elements within the quasiparticle random phase approximation (QRPA). 
In particular, it has been shown that the exact calculations for the 2v amplitude, when 
evaluated with a zero-range force, can be nicely fitted by a (1, 1)-Pade approximant of the 
form 

A^2. = ^2.(t = 0)i^, (1) 

where t = vf^/ v"^""" is the ratio between the spin-triplet strength in the particle-particle chan- 
nel and the spin-singlet strength in the pairing channel, and to and ti denote, respectively, 
the zero and the pole of M.2v and are treated as free parameters. In the same work it has 
been suggested that this result is of general validity and that any alteration of the nuclear 
hamiltonian or the configuration space cannot lead to a different functional dependence. [] 
Thus, in order to modify the theoretical predictions for the (3(3 moments in a qualitatively 
way, the QRPA itself has to be altered. The simplest amendment is to set up a particle 
number projection. This has been done by Faessler and collaborators in several recent pa- 
pers [0. They have also discussed the outcomes of the PRPA on both simple and double 
beta decays. We present here a PRPA formalism that differs in several aspects from the one 
reported by them. 

To evade the disadvantages inherent in the non-conservation of particle number, besides 
the usual proton quasiparticle transformation 

= {upCp — VpC^^ ; a|, = {upC^^ — VpC^ , 

it is useful to introduce the following canonical transformation Q 

dp = ^/Op (wpCp - Wp2;pcl) ; d* = ^ (wpcj - VpZ^c^ ; = ul + zlv^, (2) 

and likewise for neutrons. Here (Cp), with p = p, mp, p = Up £p jp and rup = rrij^, are the 

single particle creation (annihilation) operators for protons, and Cp = {—iy~^"^pcp^_mp- The 

^When the renormalization coupling constant gPP is used Q an analogous expression to the eq. (jl|) is 
vahd (with gPP's for t's). 

^ Note that the operation symbolized by star (★) does not affect the complex variable Zp. 
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vacuums for the proton operators dp is the generating function ^ 

\z^) = z^'^\BCS- p) = n (up + Zpi^pcjc^) 

p>0 

where Z is the particle number operator for protons, and \p) and \BCS;p) represent, respec- 
tively, the particle and the BCS vacuum for protons. The Bogoljubov transformation (0) 
has been introduced in the past by Ottaviani and Savoia to perform a number projection 
in spaces of two or more BCS quasiparticles. Later, their method has been used by several 
authors (e.g., ref. 0). 

First, we derive the following relations for the matrix elements of the model hamiltonian 

H: 

{BCS\HPo\BCS) = T{z\H\z), (3) 
{BCS\[A{pnJ)Hp^A^{p'n'J)f\BCS) = T {z\[B^{pnJ)H-D'^{p'n' J)]'^\z) , (4) 

where \BCS) = \BCS,p)\BCS,n) and 1^;) = |2;p)|2;n). The projection operators are: Pq = 
PzPni P+1 = -Pz+i-P/v-i, and P_i = Pz-iPn+i, with fV being the operator for projecting 
states of good particle number A/" . The transformation 




is a contour integration including the origin, and 

At(pnJ) = [ajat]^; B;{pnJ) = ya^4i-^)/24i+^)/2[d;d;;]^ 

are the proton-neutron excitation operators in the spaces of the operators a and d, respec- 
tively. 

For the norms, we get 

{BCS\Po\BCS)=No = rHo, (5) 

and 

J-'{BCS\[A{pn; J)P,A\pnJ)f\BC S) = N^{pn) = TN^ipn), (6) 
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where 

The one-body charge-exchange operators 

O^AJ) = E(P|0(^)l")4cn; O-i(J) = E(n|0(^)IP)cIcp, 

pn np 

when expressed in terms of the d operators, and after neglecting the scattering terms that 
never contribute within a two quasiparticle space, read 

0,{J) = J2[Al{pnJ)B;ipnJ) + A'_^ipnJ)B.,{pnJ)], 

pn 

with 

lO/^^ T\ _ _ T-i/^iin/' 7Mi^\ J ""p""" for /i = 1 , 



UnVp for fx = —1 . 



Al{pnJ) = -J~'{p\\OiJ)\\n) 

For their matrix elements we get 

{BCS\[A{pnJ)P^O,,{J)fPo\BCS) = 0^{pnJ)6,,,^, = T0^{pnJ)6^,^,, (7) 

with 

0,{pnJ) = {z\[B,{pnJ)0,{J)f\z) = M,{pn)Al{pnJ). 

Relations (^, (^, (|^) and reveal that, at the level of the Tamm-Dancoff approximation 
(TDA), there is a one to one correspondence between the energy spectrum of H within 
the projected basis P^A'^ {pnJ)\BC S) (relative to the ground state Po\BCS), with energy 
Eq = {BCS\hPo\BCS) /Nq), and that within the basis D*^{pnJ)\z) (relative to the ground 
state l^;), with energy Eq = {z\H\z) /Nq). As it is seen from (|^) an analogous correspondence 
also stands for the transition matrix elements 0^{pnJ) and Ofj,{pnJ). Thus within the TDA 
our method consists in the following: after evaluating the matrix elements of H and 0±i{J) 
within the d-space (< D >), the corresponding projected matrix elements (< P >) are 
obtained from the transformation: < P >= T < D >. 
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Next, the above method is extended to the RPA by following the linear response procedure 
developed by Lane and Martorell We introduce the Slater determinant 

=11(1 + x,ipnJ)Ti;ipnJ))\z), (8) 

pn-J/j, 

and if the quantities x ci-re assumed to be small, then, to second order 

= No + E [\x,{pnJ)f + \x-,{pnJ)f N-.ipn)] . (9) 

pnj/i 

The response of the system to an external time-dependent field 

is derived from the time dependent variational principle, which for the optimum Slater 
determinant solution of the perturbed problem gives 

d_ 



{6<P{z)\H - Eo + - ih-Mz)) = 0. 



If Xfj.{p^J) cire varied, keeping only the linear terms and assuming {z\0^\z) = 0, one gets 



A^X^ + Bx% + O^e-'^' - ihx^H^ = 
where the submatrices are given by 



\{pn,p'n'J) = J'\z\[B,{pnJ){H -Eol)B;{p'n'J)]''\z) 
B{pn,p'nJ) = J-\z\[D^^{pnJ)-D^{p'nJ)H]'^\z). 

Attempting a solution of the form 

X, = X,e''-'; x~, = Y;e'-\ (10) 

and performing the transformation {A^, B, Nfj_, O^) = T(A^, B, N^, 0^), we obtain 

( A,-hu;N,l B \IXA=_I OA , . 

\ B* A*_^ + hu;N*_^l ) \ yJ [Olj- ^ > 



The corresponding RPA eigenvalue problem is: 




with the eigenvectors normalized such as to obey the orthogonal relation: 

Xl{a)N,X,{P) - Y^{a)N^,Y,{P) = 
while the closure relations are: 



6af3 for E^{a) > , 
-6ai3 for Ef,{a) < 



Y: [X;{a)Y,{a)-Y.,{a)xUa)]Nl^ = 0. 

We interpret E^{aJ) as the PRPA approximation for the exact excitation energies from 
the ground state in the nucleus (Z, A^) to the states \aJ]ii) in the neighboring odd-odd 
nuclei: Efj,{aJ) = {a J; fi\H\aJ; ^)prpa — {0\H\0)prpa- Note that in contrast with the usual 
charge- conserving RPA problem, the eigenvalues do not occur in pairs ±i?^(aJ). Besides, 
both positive and the negative solutions are physically meaningful. For fi = ±1 the positive 
solutions describe excitations in the {Z±l, A^=Fl) nuclei, while the negative energy solutions 
represent the de-excitations in the (Z ± 1, A^ =i= 1) nuclei, and the positive energy excitations 
in the [Z ^ 1,N ±1) nuclei as well. Thus only one RPA equation has to be solved for 
the evaluation of single (3~ and f3~^ transitions. This is a well known feature of the charge- 



exchange modes P, 0, 10 



We can express the solution of eq. (|Tl|) in terms of the eigenvectors of (|12|) 



and hence from (^), (^) and ( PH] ) we derive the response function ||^, 

l T{<l>iz)\d,\<l>iz)) |A>J)P 



O 




2 r{<P{zmz)) tr " ''E,{aJ)-hu' 
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where 

A^(aJ) = N^'/' Y.K^PnJ)N^{pn)X;{pn- aJ) + A° ^(pnJ)iV_^(pn)y;(pn; aJ)], 

is the PRPA approximation to the matrix element of Of^{J) between the exact counterparts 
of \aJ;ii) and the ground state, i.e., A^(q;J) = {o(J;ij>\\Oh{J)\\0)prpa- 
Introducing the normahzed amphtudes, 

Xf,{pn; a J) = Nl^^{pn)Xi^{pn; a J); y^j,{pn] a J) = N]!^{pn)Y^{pn] a J), 

the PRPA equation takes the standard form 

/ A,{J) B{J) \ ( X,{aJ) X,{aJ) \ 



with 



A^,{pn,pn;J)- 



B{pn,p'n; J) 



' ^N^{pn)N^{p'n'y 

B{pn,p'n'; J) 
^N,{pn)N_,{p'n'y 

The orthogonahty and closure relations are, respectively, 
and: 



Sai3 for Ef,{a) > , 
-Saf3 for Eij_{a) < . 



E [X;{a)y,{a) - y.,{a)Xl,{a)] = 0, 

while the transition amplitudes read 

A,{aJ) = Ar-V2 J2^A'^^j^J)N'J^pn)X;{pn; aJ) + k%{pnJ)N'J^{pn)yi{wi- aJ)], (13) 

jm 
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and for J'^ = and 1^ satisfy the sum-rule [|| 

E lA.M)!'- E \A-,{aJ)\' = {2J + l){N-Z). 

The exphcit results for the projected RPA matrices are: 

A,{pn,p'n'; J) = e^-'^^'''-'-^ipn)S,^^,,^, + R',^'^^^''-'-^{pnp'n'- J), 
B{pn,p'n'; J) = R^^^ {pn.p'n'] J), 



where 



e'^^'ipn) = [R'^ip) + R',,ipp)]l'' (n) + + 

+ R'/{pn) + R\t\pn) - E^I^p)!'' {n), (14) 



are the unperturbed energies, and 

dZ tdN tdZ,N 
_ Kq Kq Kq 

^0 - + JJT + JzJn- ^^^> 
is the ground state energy. The norms are: 

No = I'l""; N,{pn) = /^-i+^(p)/^-i-^(n), 

and the remaining quantities are defined as follows: 

Roipi--) = EiW-'(m--) 

p 

+ \T.hJA'^vyAPP')l'''\pp'Pi ■ ■) + UpVpUp>Vp>g{pp')l''-\pp'pi ■ ■)], 
pp' 

Rn{PiP2 ■ ■) = ep,[ulj\pip2 ■ ■) - vlj''-\pip2 ■ ■)] 

+ j;'T.JpHKppi)Ki'-\ppm ■ ■) - vll'^-\pp, ■ ■)] 

p 

- UpVpUp^Vp^g{ppi)I^'^{ppiP2 ■ •)}, 
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Ro'^'iPi ■ -,^1 ■ ■) = j:ynvyjipn)l''''ippi ■ ■)l''-\nn, ■ ■)) 

pn 

Ri'i'ipn) = 3;\ull\pp) - vll'-\pp)]Y.3n'vl,i{pn')I 



t\ Tk'-2f >\ 

^ ' inn 



+ VKl^'inn) - vll''-\nn)]Y^j,,vl,finpV'-'ipp'), 



R22 {pn,p'n'; J) = [upVnUpiVn'I^{pp')I^'^'^{nn') + VpUnVp'Un'I^~''^{pp')I^' {nn )]Y{pn,p'n'] J) 
+ [upUnUpmn'l'^{pp')I^' {nn') + VpVnVp'Vn'l''~'^{pp')l''''^{nn')]G{pn,p'n'; J) 
Rlf{pn,p'n'] J) = I^'^{pp')I^~'^{nn') [{vpUnUp>Vn' + UpVnVp'Un')F{pn, p'n' ; J) 

- {UpUnVpiVn' + VpVnUpiUn')G{pn, p'n'] J)] . 

where Cp are the proton single particle energies (s.p.e.), i{pp') = F{pp,p'p'; J = 0), g{pp') = 
G{pp,p'p'; J = 0), f{np) = ¥{nn,pp; J = 0) and 

1 r dzy. 



2ni J z!^+^ 



The starting point of a self-consistent numerical work is the minimization of Eq with 
regard to the parameters u and v. This leads to a set of coupled equations 0, B]: 



lepUpVp - l\p[u\-vl) = 
where the quantities and Ap are defined as follows: 

= -\jp^Y.h'Up'Vp'g{pp)I^~\p), 
p' 

Bp = epl'-\p) + E vl,i{pp')l'-\pp') 

p' 



(16) 



I -'- "■-2 "2 2 - 



Vp{p')D'-\pp') - VpD'-\p') 



Z-2i 



1^ 



Vp{plf)D^-\pp^f) - VpD^-\p)' P > 



1^ 



+ Up,Vp,UpnVpng{p)f) 



Vp{p'p")D'-\pp'p")-UpD'-\py 



1^ 
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with 

n 

D\p..) = l\p..)-I'^^\p..). 
To inquire into the effect of the number projection method we have done numerical 



calculations for several 2z//3/3 emitting nuclei, following the procedure outlined in ref. |T2 |. 
This implies the employment of the s.p.e. from the neighboring odd-odd nuclei, and therefore 
all self-energy corrections in eqs. ([T^), (pIS]) and (p!6D, for both like and unlike particles, have 
been neglected. Besides, and for the sake of comparison with results presented in ref. 0, 
the 0^ 0+ 2z/ transition amplitude has been approximated as 

-^2. = 77 E 



2 ^, \E^=,{aJ = 1) + A E^=_i(a'J = 1) - A, 



xA^=i(«J = l){aJ =l;fi = l\a'J = = -l)A^=_i(a'J = 1), (17) 

where unbarred and barred quantities indicate that the quasiparticles and excitations are 
defined with respect to the initial (Z, N) nucleus, and final {Z+2, N—2) nucleus, respectively, 
and A = (Eo - ^o) /2. 

As an example, the PQRPA and QRPA calculations for the ^^Ge ^^5*6 transition 
are displayed in fig. [l|. Also, in the same figure it is shown the result of a hybrid model 
calculation, in which the PQRPA is used, but with the BCS parameters u and v employed 
within the QRPA. It can be observed that all the three results are qualitative similar, in the 
sense that the f3(3 moments always vary rather abruptly in the physically relevant interval 
to^t^ti, and are well represented by eq. (|^); the corresponding parameters to and ti that fit 
the exact results, together with M.2v{t = 0), are listed in table 0. The single mode model for 



the (5(3 decay [||, suggests that the dominant effect of the projection method is to insert, 
into the numerator and the denominator of eq. (|lD, terms of the order of ^ + ^) 5 where 
2Vtp and 2f2„ are, respectively, the degeneracy of the proton and neutron levels. 

We end the present letter by mentioning the main differences between the PRPA formal- 
ism reported in ref. p[ and the one presented here: 



9 



a) They solve two different RPA equations for tlie {Z + 1, N — 1) and {Z — 1, + 1) nuclei. 
Contrarily, we only deal with one RPA problem, in which the l3~^ spectrum is viewed as an 
extension of the P~ spectrum to negative energies ||, ^, [10[ . 

b) Their matrices A and B are defined through the method of commutation relations. Yet, 
we feel that this technique cannot be applied with the projected procedure, as it automati- 
cally excludes the _Ro terms in the expression for the unperturbed energies. But these 
terms in no way cancel among themselves P], and thus, our unperturbed energies are 
necessarily distinct from those obtained in ref. 0. (It is worth noting that the commutation 
method has neither been used in the derivation of the PRPA for like particles |1TT| .) 



c) Their result for the matrix B also disagrees with our result and we do not find any expla- 
nation for such a discrepancy. In our formalism the ground state correlations always occur 
in the {Z, N) nucleus. 

d) We do not have to deal with spurious components in the wave functions of the intermedi- 
ate double-odd nucleus as it apparently happens in ref. 0, and therefore A^=i, A^=_i and 
B are square matrices of dimension equal to the number of proton-neutron states with a 
given spin and parity. 

e) The matrix element for the one-body transition operator, given by eq. {^^, has not been 
explicitly shown in their previous works. 

In summary, the QRPA and PQRPA yield qualitatively similar results, and, in spite of 
above mentioned differences between the formalism employed by the Tiibingen group and 
ours, the numerical results that we get for the moments A^2i/ are not in essence distinct from 
those obtained in ref. 101. 
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Figure Captions 



Figure 1: Calculated double beta decay matrix elements A42u (in units of [Me\^]^^) for 
^®Ge, as a function of the particle-particle S = 1, T = couphng constant t. Solid and dotted 
curves correspond to the projected (PQRPA) and unprojected (QRPA) results, respectively. 
For comparison a hybrid model result, in which the PQRPA is used with parameters u and 
V determined as in the QRPA, is also presented (dashed curve). 



Tables 

Table 1: The coefficients A42u{t = 0), to and ti for ''^Ge in the parametrization of 2z/ moments 
by eq. (p. The exact curves, shown in fig. are not distinguished visually from the fitted 
ones with the parameters listed here. 





QRPA 


PQRPA 


PQRPA 
(hybrid) 


-M2,it = 0) 


0.385 


0.370 


0.297 


to 


1.170 


1.038 


0.950 
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1.591 


1.447 


1.421 
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